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Bearing-Based Formation Stabilization with Directed Interaction Topologies 

Shiyu Zhao and Daniel Zelazo 


Abstract —This paper studies the problem of stabilizing 
target formations specified by inter-neighbor bearings with 
relative position measurements. While the undirected case has 
been studied in the existing works, this paper focuses on the 
case where the interaction topology is directed. It is shown that 
a linear distributed control law, which was proposed previously 
for undirected cases, can still be applied to the directed case. 
The formation stability in the directed case, however, relies 
on a new notion termed bearing persistence, which describes 
whether or not the directed underlying graph is persistent with 
the bearing rigidity of a formation. If a target formation is 
not bearing persistent, undesired equilibriums will appear and 
global formation stability cannot be guaranteed. The notion of 
bearing persistence is defined by the bearing Laplacian matrix 
and illustrated by simulation examples. 

I. Introduction 

In recent years, there has been a growing research interest 
on bearing-based formation control [ 1 ]—[8] and bearing- 
based network localization [9]—[13]. The bearing-based ap¬ 
proaches provide solutions to control or localize networks in 
arbitrary dimensional spaces merely using the inter-neighbor 
bearing information. As shown in [7], the bearing-based 
approach also provides a simple solution to formation scale 
control, which is complicated to solve using the distance or 
relative position based approaches [14], [15]. While most of 
the existing works focused on the undirected case, this paper 
focuses on bearing-based formation control with directed 
interaction topologies. 

The specific problem considered in this paper is described 
below. Suppose the target formation is static and its desired 
geometric pattern is defined by the inter-neighbor bearing 
constraints. Each agent can obtain the relative positions of its 
neighbors via wireless communication. The control objective 
is to steer the formation from an initial configuration to a 
final configuration that satisfies the inter-neighbor bearing 
constraints. The linear formation control law considered in 
this paper is the one proposed in our previous work [6] for 
formations with undirected interaction topologies. It will be 
shown in this paper that the control law can also be applied 
to the directed case, but the formation stability in the directed 
case is more complicated to analyze. 

In the undirected case, the formation stability merely relies 
on the bearing rigidity properties of the target formation [6]. 
The bearing rigidity of networks or formations has been 
studied in [2], [3], [5], [8]. It is shown in [5] that the 
geometric pattern of a target formation can be uniquely 
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determined by the inter-neighbor bearings if and only if 
the target formation is infinitesimally bearing rigid. It is 
further shown in [6] that the global formation stability will be 
guaranteed if the target formation is infinitesimally bearing 
rigid. However, in the directed case, the formation stability 
will rely on not only the bearing rigidity properties but also 
a new notion termed bearing persistence as observed in this 
paper. The notion of bearing persistence is defined based on 
a special matrix, termed the bearing Laplacian , which plays 
important roles in the bearing-based formation control and 
network localization problems as shown in [6], [7], [13]. 
The bearing Laplacian of a formation can be viewed as a 
matrix-weighted graph Laplacian. In the undirected case, the 
bearing Laplacian is symmetric positive semi-definite and 
its null space is exactly the same as that of the bearing 
rigidity matrix. However, the bearing Laplacian may lose 
these elegant properties in the directed case. 

The work presented in this paper is a first step towards 
solving the general problem of bearing-based formation 
control with directed interaction topologies. The notion of 
bearing persistence is defined and explored. Specifically, a 
formation is bearing persistent when the null spaces of the 
bearing Laplacian and the bearing rigidity matrix are the 
same. We further show that if the formation is not bearing 
persistent, undesired equilibriums will appear and the global 
formation stability cannot be guaranteed. We also present 
a sufficient condition to ensure bearing persistence for for¬ 
mations in the two-dimensional space. The spectrum of the 
bearing Laplacian in the directed case, which determines the 
formation stability, is however not analyzed in this paper and 
will be addressed in the future. Linally, since the bearing- 
based formation control problem is mathematically equiv¬ 
alent to the bearing-only network localization problem as 
shown in [6], the results presented in this paper can be easily 
transferred to solve the bearing-only network localization 
problem with directed interaction topologies. 

Notations: Given A t e R pX! for i = 1 ,n, denote 
diag(Aj) = blkdiag{Ai, ... ,A n } e R” pX119 . Let Null(-) 
and Range(-) be the null space and range space of a matrix, 
respectively. Denote I,i £ R. dxd as the identity matrix, and 
1 = [1,..., 1] T . Let || • || be the Euclidian norm of a vector 
or the spectral norm of a matrix, and <8> be the Kronecker 
product. 

II. Preliminaries to Bearing Rigidity Theory 

Bearing rigidity theory plays a key role in the analysis of 
bearing-based distributed control and estimation problems. 
In this section, we revisit some notions and results in the 
bearing rigidity theory [5], It is worth noting that the bearing 



rigidity theory presented in [5] is developed for the case 
of undirected graphs. However, it can be generalized to the 
directed case without any difficulties. 

We first introduce an orthogonal projection matrix operator 
that will be used throughout the paper. For any nonzero 
vector x € R d (d > 2), define P : —> R dxd as 



For notational simplicity, we denote P x = P{x). The matrix 
P x geometrically projects any vector onto the orthogonal 
compliment of x. It is easily verified that Pj = P x , P;. = 
P x , and P x is positive semi-definite. Moreover, Null ( P x ) = 
spanja;} and the eigenvalues of P x are {0, Any two 

nonzero vectors x, y G R d are parallel if and only if P x y = 0 
(or equivalently P y x = 0). 

Given a finite collection of n points {p,:}” =1 in R. d (n > 2, 
d > 2), denote p = [pj,... ,p^] T G K dn . A formation in 
R d , denoted as G(p ), is a directed graph Q = (V, £) together 
with p, where vertex i € V in the graph is mapped to the 
point pi. For a formation G(p), define the edge vector and 
the bearing , respectively, as 

e ij ~ Pj -Pi> 9ij ~ e ij/|| e ij||% V(«, j) G £. 

The bearing g tJ is a unit vector. 

Definition 1 (Bearing Equivalence). Two formations G (p) 
and G(p') are bearing equivalent if P gij g£ = 0 for all 

ihj) e £■ 

By Definition 1, bearing equivalent formations have par¬ 
allel inter-neighbor bearings. Suppose |£| = to and index all 
the directed edges from 1 to to. Re-express the edge vector 
and the bearing as eu and gk = efc/||efc||, V/c G {1,... ,to}. 
Let e = [ej 1 ..., e^] T and g = [gj,..., g^] T . Note e 
satisfies e = Up where H = 11 0 1,i and H is the incidence 
matrix of the graph Q [16]. Define the bearing function 
F b : R dn -> R dm as 

Fb(p) — [gj 5mJ T - 

The bearing function describes all the bearings in the forma¬ 
tion. The bearing rigidity matrix is defined as the Jacobian 
of the bearing function, 

R B (p) = dF ^ G (1) 

op 

Let Sp be a variation of p. If J?b(p)^P = 0, then 5p is called 
an infinitesimal bearing motion of G(p). Any formation 
always has two kinds of trivial infinitesimal bearing motions: 
translation and scaling of the entire formation. We next define 
one of the most important concepts in bearing rigidity theory. 

Definition 2 (Infinitesimal Bearing Rigidity). A formation 
is infinitesimally bearing rigid if all the infinitesimal bearing 
motions of the formation are trivial. 

Some useful results are listed as below; for proofs, please 
refer to [5]. 


Lemma 1 ([5]). For any formation G(p), the bearing rigidity 
matrix defined in (1) satisfies 

(a) R B {p) = dmg{P gk /\\e k \\) H; 

(b) rank(J? B ) < dn—d—1 and spanjlgi/^p} C Null(J? B ). 

Theorem 1 ([5]). Any two formations G(p) and G{p') are 
bearing equivalent if and only if Rb(p)p' = 0. 

Theorem 2 ([5]). For any formation G ( p ), the following 
statements are equivalent: 

(a) G(p) is infinitesimally bearing rigid; 

(b) Gip) can be uniquely determined up to a translation and 
a scaling factor by the inter-neighbor bearings; 

(c) rank (Rb) = dn — d — 1; 

(d) Nu 11(.R b ) = spanjl ® Id,p}- 

III. Problem Formulation 

In this section, we first state the problem of bearing-based 
formation stabilization, and then present a linear bearing- 
based formation control law. 

A. Problem Statement 

Consider a formation of n agents in R. d (n > 2, d > 2). 
Denote p, G R d as the position of agent i G {l,...,n}, 
and let p = [pj, ■ ■ ■, Pn] T G R dn . The dynamics of each 
agent is a single integrator, Pi(f) = ufit), where ufit) G 
W 1 is the input to be designed. The underlying graph Q = 
(V, £) is directed, fixed, and connected. If (i.j) G £, we 
say agent i can “see” agent j, which means agent i can 
access the relative position of agent j. We call G(p) as a 
directed (respectively, undirected) formation if G is directed 
(respectively, undirected). 

The aim of the control is to stabilize a target formation 
specified by constant bearing constraints {(]*;}(;,j)e£- The 
bearing constraints must be feasible such that there exist 
formations satisfying these constraints. The problem to be 
solved in this paper is stated as below. 

Problem 1 . Given feasible constant bearing constraints 
an d an initial configuration p(0), design ufit) 
(t G V) based on relative position measurements {p t (t) — 
Pj(t)}jeAfi suc h ii mt Pg* 9ij(t) —> 0 as t —> oo for all 

ihj) e £■ 

Let Gip*) be an arbitrary formation that satisfies the 
bearing constraints {g*j}(ij)e£- Problem 1 actually requires 
that the formation G{p{t)) converges to a final formation 
that is bearing equivalent to G{p*)- This means that we 
allow the bearings g., 3 (t) to converge to either gf or — gT. 
This is different from the problem setup in [5] where gij(t) 
must converge to gf. An illustrative example is shown in 
Figure 1. The formation in Figure 1(a) has the bearings as 
[g*j}(i,j)e£, whereas the formation in Figure 1(b) has the 
opposite bearings as {— 9ij}(i,j)e£- Problem 1 allows a final 
formation to be either (a) or (b). 





Fig. 2: The geometric interpretation of control law (2). 


B. Bearing-Based Formation Control Law 

The bearing-based formation control law considered in this 
paper is 

Pi{t) = ~Y1 P 9h ~ (*)) ’ * e V ’ (2) 

je/Vj 

where P g *. = Id, — {g*j){9ij) T - Control law (2) was originally 
proposed in [6] for formations with undirected interaction 
topologies. As will be shown in this paper later, it can also be 
applied to the directed case. Several remarks on the control 
law are given below. Firstly, control law (2) is distributed 
because the control of each agent only requires the relative 
positions of its neighbors. Secondly, the matrix P g *. is an 
orthogonal projection matrix and hence control law (2) has a 
clear geometric meaning. As shown in Figure 2, the control 
term — P g *.(pi(t ) — Pj(t)) is the orthogonal projection of 
Pj(t) — Pi(t) onto the orthogonal compliment of and 
hence it acts to reduce the bearing error. 

The matrix expression of control law (2) is 

P(t) = ~ L B p(t ), (3) 

where Lb G R dnxdr! and the ijlh block submatrix of Lb is 

[Cb]*? = Odxdj i 7^ jt (hj) ^ 

[-^-'B] = ~ p> g*ji * 7^ ji (*> j) C £, 

[Lr]u = P 9*j’ * e ^ 

The matrix Lb can be interpreted as a matrix-weighted graph 
Laplacian. It is jointly determined by the topological and 
Euclidean structure of the formation. We call Lb the bearing 
Laplacian since it carries the information of both the bearings 
and the underlying graph of the formation. The spectrum and 
null space of Lb fully determine the stability of the formation 
dynamics (3). 

IV. Bearing Persistence and Formation Stability 

In this section, we define the notion of bearing persistence 
by the bearing Laplacian and analyze the impact of the 
bearing persistence on the global formation stability of 
system (3). 


We note that three subspaces play key roles in the forma¬ 
tion stability analysis: NuII(Lb), Nu11(.Rb), and span{l ® 
Id,p*}- First, Null (Lb) is the space where the solution p(t) 
converges when the eigenvalues of L B are in the closed right 
half of the complex plane. Second, Nu11(I?b) is the space 
of all the formations that are bearing equivalent to the target 
formation according to Theorem 1. The equality NuII(Lb) = 
NuII(^b) is required in order to solve Problem 1. Third, 
spanjl ® I<uP*} is the space of all the formations that 
have the same shape as the target formation (but can have 
different translations and scales). The equality NuII(Lb) = 
spanjl ® Id,P*} is desired when the formation is required 
to converge to a formation that has the same shape as the 
target formation. We next explore the relationship between 
the three subspaces. 

We first consider the simple and special case where the 
graph is undirected. An undirected graph can be viewed as 
a special directed graph where (i,j) £ f (j, i) £ 8. In 
the undirected case, the properties of the bearing Laplacian 
can be easily obtained as below. The proof can be found in 
[13], For the sake of completeness, we also provide a proof 
below. 

Theorem 3. For an undirected formation Q (p), the bearing 
Laplacian Lb is symmetric positive semi-definite and satisfies 

NuU(Lb) = Null(-Rfl). 

Proof Each undirected edge in the graph can be viewed 
as two directed edges with opposite directions. Index the 
directed edges from 1 to m. Then Lb can be written as 
Lb = H T dia,g(P gk )H, which is clearly symmetric and 
positive semi-definite. Moreover, due to P gk = Pj k P gk , L B 
can be rewritten as 

L b = g T diag(P g T fc ) diag(P 3fc )Lf . 

kl Rb 

Note Rb = diag(11 e*; 11) L?b where Rb is the bearing rigidity 
matrix. As a result, NullfL’/j) = NuII(L’b) = NuII(Lb). □ 

Based on Theorem 3, the solution p(t) to system (3) 
in the undirected case converges into NuII(L’b) and hence 
the finally converged formation is bearing equivalent to the 
target formation Qfp*). In order to make p(t) converge 
to spanjl ® Id,p *}, it is required the formation to be 
infinitesimally bearing rigid according to Theorem 2. 

We now consider the general directed case. When the 
graph is directed, the bearing Laplacian is not symmetric any 
more and we do not have NuII(Lb) = NuII(L’b) in general. 

Theorem 4. For a directed formation Qip), the bearing 
Laplacian Lb satisfies 

spanjl ® Id,p} Q Null(-Rfl) C Null(L B ). 

Proof The fact spanjl ® Id,p} Q Nu11(I?b) is given in 
Lemma 1(b). 

We next prove Nu11(.Rb) f NuII(Lb). In order to do that, 
we firstly show that for any vector p> £ the equality 




Rbp' = 0 holds if and only if 

Pg i M-ti)= °> (4) 

As shown in Lemma 1, by indexing the directed edges 
from 1 to m, the bearing rigidity matrix can be expressed 
by R b = dia.g(P gk /\\e k \\)H. As a result, R B p' = 0 <G> 
diag(P 9fc /11e*;11) Hp' = diag(P^ fc /||e fe ||)e' = 0 P 9k e' k = 
0, Vfc G {1,... ,m}. Therefore, R B p' = 0 if and only if (4) 
holds. We secondly show that the equality L B p' = 0 holds 
if and only if 

Y J p ai M~p' j ) = 0, Vi G V. (5) 

jeK 

Since 

LbP = YjeUi Pgij (Pi ~ Pj ) , 

we have L B p' = 0 if and only if (5) holds. Finally, since 
equation (4) implies equation (5), any vector in NuII(Pb) is 
also in Null (Lb) and hence NuII(Pb) C NuII(Lb). □ 

If NuII(Lb) is larger than NuII(L’b), the formation may 
converge to a final formation in Null (Lb) but not in 
NuII(Pb), which is not bearing equivalent to the target 
formation. Hence it is required that NuII(Pb) = NuII(Lb) 
in order to solve Problem 1, which motivates the following 
notion. 

Definition 3 (Bearing Persistence). A directed formation 
Q{p) is bearing persistent if Null (Lb) = Null(L#). 

The name of persistence actually is inherited from the 
notion of persistence in distance-based formation control 
problems [17], [18]. Bearing persistence describes whether 
or not the directed graph is persistent with the bearing 
rigidity of a formation. If they are not persistent, undesired 
equilibriums would appear for system (3) and the target 
formation cannot be globally stabilized. In addition, it fol¬ 
lows from Theorem 3 that undirected formations are always 
bearing persistent. The next result immediately follows from 
Theorem 4 and Definition 3. 

Theorem 5. For a directed formation Q{p), the equality 
NuII(Lb) = span{ 1 ® I d ,p} 

holds if and only if the formation is both infinitesimally 
bearing rigid and bearing persistent. 

Proof According to Theorem 4, the equality NuII(Lb) = 
spanjl <g> I d ,p} holds if and only if NuII(Lb) = NuII(Lb) 
and NuII(-Rb) = spanjl ® I d ,p }, where the former one 
holds if and only if the formation is bearing persistent by 
definition and the latter one holds if and only if the formation 
is infinitesimally bearing rigid by Theorem 2. □ 

Theorem 5 is important in the sense that it gives the nec¬ 
essary and sufficient condition for the equality NuII(Lb) = 
spanjl® I d ,p} in the directed case. This equality is desired 



Fig. 4: Examples of bearing non-persistent formations 


when the formation is required to converge to a formation 
that has the same shape as the target formation. 

A. Illustrative Examples 

Examples are shown in Figures 3 and 4 to illustrate the 
notion of bearing persistence. We next specifically examine 
the two formations in Figure 3(a) and Figure 4(a). For the 
two formations, the agents have exactly the same positions 
as pi = [—1,1] T , p 2 = [1,1] T , P3 = [1, —1] T - and p 4 = 
[— 1, — 1] 1 - The underlying graphs of the two formations, 
however, are different; the edges between agents 1 and 2 
have the opposite directions. First of all, both of the two 
formations are infinitesimally bearing rigid because it can be 
verified that for both of the formations rank(i?B) = 2n — 3. 
The null space of R B for the two formations are the same 
and given by 


Null (Lb) = span ■ 


= span{ 1 <g> h,p}- 


For the formation in Figure 3(a), it can be verified that 
Null(L B ) = NuII(L’b). As a result, this formation is bearing 
persistent. However, for the formation in Figure 4(a), the null 
space of Lb is 


NuII(Lb) = span • 


0 

1 

1 

2 

1 

0 

1 

-1 

0 

1 

1 

1 

1 

0 

1 

-2 

0 

1 

-1 

0 

1 

0 

-1 

0 

0 

1 

-1 

0 






As can be seen, in addition to the basis vectors in Nu11(_Rb), 
an extra vector (the fourth) is contained in NuII(Lb)- As a 
result, this formation is not bearing persistent. We will later 
show by simulation that this non-persistent formation cannot 
be globally stabilized. 

It must be noted that the bearing persistence of a formation 
is independent to the bearing rigidity of the formation. For 
example, the directed formation as shown in Figure 3(a) is 
both bearing persistent and infinitesimally bearing rigid. The 
one in Figure 3(b) is bearing persistent but not bearing rigid. 
The one in Figure 4(a) is infinitesimally bearing rigid but 
not bearing persistent. The one in Figure 4(b) is neither 
infinitesimally bearing rigid nor bearing persistent. 

B. A Sufficient Condition for Bearing Persistence 

An important problem that follows the definition of bear¬ 
ing persistence is what kind of formations are bearing 
persistent. We next give a sufficient condition for bearing 
persistence. All the persistent formations given in Figure 3 
can be examined by this sufficient condition. For a directed 
edge (i,j) £ £, since this edge leaves agent i and enters 
agent j, we call it an outgoing edge for agent i and an 
incoming edge for agent j. In the formation dynamics (2), 
each agent is responsible for achieving all the bearings of 
its outgoing edges. It is intuitively reasonable that when 
the number of the outgoing edges of an agent is large, the 
chance that the agent is not able to archive all the desired 
edge bearings is high. This intuition is consistent with the 
following sufficient condition for bearing persistence in R 2 . 

Proposition 1. For a directed formation G(p) in R , if each 
agent has at most two non-collinear outgoing edges, then the 
formation is bearing persistent. 

Proof. We prove that in this case any vector x £ R 2ra in 
NuII(Lb) is also in NuII(-Rb)- Note L^x = 0 if and only if 
equation (5) holds. Suppose agent i £ V has two outgoing 
edges (i.e., two neighbors). Let the two neighbors be agents 
k and I. Then, equation (5) implies 

Pgik (*R ^- k ) T ('R tcf) 0. (6) 

Note the orthogonal projection matrices P g . k and P gu project 
vectors onto the orthogonal vectors of g,k and gn, respec¬ 
tively. Since g,t- and g,p are not collinear, equation (6) holds 
if and only if P gik {xi — Xk) = 0 and P gu (xi — xf) = 0. 
The case where agent i has merely one outgoing edge (i.e., 
one neighbor) is trivial to analyze. In summary, if each agent 
has at most two non-collinear outgoing edges, equation (5) 
implies (4) and consequently NuII(Lb) = NuII(-Rb)- □ 

The sufficient condition in Proposition 1 provides an 
intuitive way to check the bearing persistence of a formation 
in R 2 . All the formations in Figure 3 can be concluded as 
bearing persistent by Proposition 1. It should be noted that 
the outgoing edges in Proposition 1 must be non-collinear; 
otherwise, the formation may be not bearing persistent. For 
example, although in the formation given in Figure 4(a) each 
agent has at most two outgoing edges, the formation is still 


non-persistent because the two outgoing edges of agent 3 are 
collinear. 

C. Discussion of the Formation Stability 

Up to now, we have explored the notion of bearing 
persistence and showed that the global stability cannot be 
guaranteed for formations that are not bearing persistent. An 
important problem that has not been discussed is whether all 
the eigenvalues of the bearing Laplacian are in the closed 
right half of the complex plane. For undirected formations, 
the answer is positive as shown in Theorem 3. For directed 
formations, the problem is nontrivial to solve and will be 
the subject of future work. But according to a large amount 
of simulations, it is conjectured that the eigenvalues of the 
bearing Laplacian of a directed formation have nonnegative 
real parts. If the conjecture is correct, then the formation 
dynamics (3) is marginally stable. Based on the linear 
system theory, the solution p(t) converges to the projection 
of p(0) onto NuII(Lb)- If the target formation is bearing 
persistent, the final formation is in Nullf/jR ) and hence 
Problem 1 can be solved. Furthermore, if the target formation 
is both bearing persistent and infinitesimally bearing rigid 
(i.e., NuII(Lb) = span{l„ (g) Id}), the final formation has 
exactly the same shape (but maybe different translation and 
scale) as the desired target formation. 

V. Simulation Examples 

We firstly consider the bearing persistent target formation 
given in Figure 3(a). The simulation results are given in Fig¬ 
ure 5. As can be seen, the control law (2) can steer the forma¬ 
tion to a final formation that is bearing equivalent to the target 
formation. We secondly consider the target formation given 
in Figure 4(a) that is not bearing persistent. The simulation 
results are given in Figure 6. As can be seen, the control law 
(2) steers the formation to an undesired formation. That is 
because agent 2 is not able to find a position to satisfy all the 
three bearing constraints g^i, g^ 4 , and <? 23 - The control input 
||p(f)|| = ||LB_p(f)|| converges to zero, which means p{t) 
converges into NuII(Lb), but as shown in Figure 6(d) the 
bearing error J2(i,j)e£ WPg^iPiitf-Pj (f))ll can not converge 
to zero, which means p(t) does not converge to NuII(-Rb)- 
As a result, p(t) finally converges into NuII(Lb)\Nu11(.Rb)- 

VI. Conclusions and Future Work 

In this paper, we studied bearing-based formation sta¬ 
bilization with directed interaction topologies. We studied 
a linear distributed control law that can stabilize bearing- 
constrained target formations in arbitrary dimensions. It was 
shown that an important notion, bearing persistence, emerges 
and critically affects the formation stability in the directed 
case. When a target formation is not bearing persistent, 
undesired equilibriums will appear and the global formation 
stability cannot be guaranteed. A sufficient condition that 
provides an intuitive way to examine the bearing persistence 
of two-dimensional formations was also proposed. 

The work presented in this paper is a first step towards 
solving the general problem of bearing-based formation con¬ 
trol. Several important problems are not completely solved 
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Fig. 5: The simulation result for the bearing persistent target formation as shown in Figure 3(a). 
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Fig. 6: The simulation result for the bearing non-persistent target formation as shown in Figure 4(a). 



in this paper and should be studied in the future. First, 
the spectrum of the bearing Laplacian with directed graphs 
is not solved by this paper. Second, although a sufficient 
condition for the bearing persistence of formations in the 
two-dimensional space was proposed, the necessary and 
sufficient conditions have not been obtained. Third, the left 
null space of the bearing Laplacian is also important for 
determining the final converged value of the formation and 
should be studied. Finally, this paper only considered the 
leaderless scheme. In order to control the centroid and scale 
as well as the maneuver of the formation, the leader-follower 
scheme should be studied in the future. 
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